Evolution of Large Scale Curvature Fluctuations 
During the Perturbative Decay of the Inflaton 

Fabrizio Di MarccQ 

Dipartimento di Fisica, Universita degli Studi di Bologna and INFN Via Irnerio, 46-40126, Bologna, Italy 

Fabio FinellQ 

INAF/IASF-BO, Istituto di Astrofisica Spaziale e Fisica Cosmica di Bologna 
Via GobetU 101, 1-40129 Bologna - Italy 
INAF/OAB, Osservatorio Astronomico di Bologna, Via Ranzani 1, 1-40127 Bologna - Italy and 
INFN, Sezione di Bologna, Via Irnerio 46, 1-40126 Bologna, Italy 

Alessandro Gruppusc0 

' INAF/IASF-BO, Istituto di Astrofisica Spaziale e Fisica Cosmica di Bologna 

I Via Gobetti 101, 1-40129 Bologna - Italy and 

^S) ■ INFN, Sezione di Bologna, Via Irnerio 46, 1-40126 Bologna, Italy 

^ • (Dated: February 5, 2008) 

OWe study the evolution of cosmological fluctuations during and after inflation driven by a scalar 

, field coupled to a perfect fluid through a friction term. During the slow-roll regime for the scalar 

. field, the perfect fluid is also frozen and isocurvature perturbations are generated. After the end 

^—1 ' of inflation, during the decay of the inflaton, we find that a change in the observationally relevant 

I large scale curvature fluctuations is possible. 

1^ • PACS numbers: 98.80.Cq 
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^ ; I. INTRODUCTION 

, The inflationary paradigm provides an explanation for the large-scale curvature fluctuations seen in the pattern of 
J> ■ anisotropics of the cosmic microwave background and of the large scale structure [l| . In an inflationary scenario driven 
' by a single scalar field, large-scale curvature perturbations are generated by the amplification of quantum fiuctuations 
OO of the infiaton field during the accelerated expansion era and remain constant until their reentry in the Hubble radius 
^""^ in the absence of post-inflationary changes. 

In the context of cold inflation, the release of entropy is left for the post-inflationary evolution, during which the 
lO inflaton decays into other intermediate flelds or directly into the matter our present universe is made of. 

[ The possibility of a change in the amplitude of large scale curvature perturbations during preheating has been 

■ widely investigated [1, H, 0, 11,J1, 0, Such a change is possible due to the coupling of isocurvature and curvature 
perturbations on large scales [3, ■ For additional scalar flelds coupled to the inflaton, the possibility of a post- 
inflationary ampliflcation of curvature perturbations is very model dependent, but possible and strongly related to 
the spectrum of the decay products generated during inflation. 

Although a post-inflationary change of large scale curvature perturbations during reheating might be an unexpected 

■ twist for successfull single field inflationary models, it may be useful for the single held models whose parameters 
clash with observations. Along this line, after it was shown that curvature perturbations can indeed change during 
preheatin g, it was proposed that large scale curvature perturbations may be seeded by fight scalar fields during 
reheating [ll|. The proposal to seed curvature perturbations at reheating requires three components, according to 
pd| : the inflaton, the decay products (taken as a perfect fluid) coupled to the inflaton and a light modulus on which 
the decay rate depends. 

Although a lot of attention has been paid to system of scalar fields, also interacting through the kinetic terms 
[S > here instead we analyze the evolution of cosmological perturbations in an inflationary scenario driven by a 
standard scalar fleld coupled to a perfect fluid (with equation of state pp = ojfPf) through a friction term F. This 
setting is motivated from the old theory of reheating [13, [H, [3] to describe the decay of the inflaton into matter 
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after the accelerated stage, but we believe it is also interesting for several reasons. First, fluctuations of a scalar 
field with a non-vanishing potential are different from those of a perfect fluid. As a second point, it is conceivable 
to explore the consequences of a coupling for the inflaton not dictated by its interaction with another scalar field. 
Third, the interaction between the inflaton and the perfect fluid is invisible to gravity, in contrast with what happens 
to scalar fields where a term g^cfP'x^ appears as part of the potential driving the inflationary stage. For these three 
reasons duplication of the results of two-field inflationary models is not expected. This formalism is also suitable for 
investigating the prediction of warm inflation [13, [isl . , in which the dissipative term F is much larger than the 
Hubble parameter during inflation. Note that here the numerical value for the ratio T / H is left free in this paper. 

We stress that the main goal of this paper is to study cosmological perturbations for an inflaton coupled to a perfect 
fluid during and after inflation. To our knowledge, in previous studies of cosmological perturbations in warm inflation 
the post-inflationary stage was not considered [l^, [la]- Apart from warm inflation, this setting has been applied 
only after inflation to describe the reheating stage with a decay rate computed by quantum field theory methods. 
However, the inflationary geometric amplification of zero-point fluctuations applies to any component present in the 
Lagrangian. If we think about the inflaton coupling to a perfect fluid as an effective way to describe its couplings to 
other degrees of freedom present in the Lagrangian, it is conceivable to consider the perfect fluid contribution also 
during inflation. 

The outline of paper is as follow. In section II and III we present the governing equations for the background 
and scalar perturbations in the longitudinal gauge, respectively. In section IV we present the evolution of scalar 
perturbations in the uniform curvature gauge (UCG henceforth) , from which is easy to obtain the coupled equations 
for the Mukhanov variable (the gauge invariant scalar field fluctuation [13]) and the Lukash variable (the gauge 
invariant fluid fluctuation [21!]). The numerical analysis presented is based on this last set of equations. In section V 
we conclude. 



II. BACKGROUND EVOLUTION 



In this section we present the background equations for the scalar field plus perfect fluid and several background 
quantities, as the speed of sound for the various components. The equations of motion for the scalar field and the 
fiuid we consider are the following: 



-?>H<t) -T4,-V4, 



(1) 



(2) 



where H = a/ a, wj? is the state parameter of the fluid - which we consider non negative - F is the friction coefficient - 
which can depend on any other quantity. Note that this coupling has a fully covariant description [19] and that other 
coupling of the inflaton to a perfect fluid may be considered 11[. The above equations correspond to: 



where i = (j), F and 



-X, 



Pi + 3H{1 + Wi)pi = Xi , 
-TcfP. The Einstein equations are: 



+ ¥{(/)) + PF 



8nG 



-Ptot ■ 



(3) 



(4) 



H 



-4nG 



0^ + Pf(1 + (^f) = -47rG' (ptot +Ptot) = -47rGptot (1 + Wtot) , 



(5) 



where it is clear that the dissipative term F does not enter explicitly in the Einstein equations. The dissipative term 
F slows down the inflaton in a slow-roll regime and damps it into the fluid during the oscillatory stage when the value 
of the Hubble rate H drops below F. 

For future convenience we define the parameter uj^ as: 



P4, 4>y2 - V{<P) 
= — = 



P0 (/.V2 + v{4>) 



and the sound of speed associated to scalar field as: 



2V^ 



3iJ0(l-t-r) 



(6) 



(7) 
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FIG. 1: Typical evolution of background quantities from inflation to radiation (wp = 1/3 , F/m — 0.5) for V{(f>) = m^fj)^ /2 as 
function of an adimensional cosmic time mt. On the left panel we test validity of Eq. (|10p : we show the evolution of </>\/G 
in the coupled (solid line) and uncoupled (dashed line) case. Note how Eq. (|10|) is a good approximation during inflation. In 
the second panel the time evolution of / {iHrn?) (solid line) and 47rGpF(l + WF)/m^ (dashed line). The time evolution 

of p^l ptot (solid) and pp / ptot (dashed) is displayed in the third panel. On the right panel it is shown the evolution of —H/H^ 
(solid line) and its contribution from the scalar field 4TTG(j)^ / (long-dashed line) and from radiation 47rGpF(l + wp)/H'^ 
(short dashed line). 



where r = r/(3if). This last formula is the generalisation of Eq. (6) of [111 to the case r 7^ 0. We also define the 
total speed of sound: 



■-tot 



P4>+PF Wtot 2 , PF + PF - r{p^ + P4,) , 2\ ^QN 



P4> + PF 3iJ(l + Wtot) Ptot + Ptot 

When the inflaton slow-rolls, the following approximations hold 



3H{1 + r) 



Pf 



■iH{l + ujF) {1 + ^f) 



(10) 



It is clear that the perfect fluid also remains almost frozen because of the coupling to the inflaton, i.e. pp ~ 0. The 
(cosmic) time derivative of the Hubble distance is: 



H 



(11) 



where ~ denotes the validity during slow-roll. 

We note that this coupling of the inflaton </) to a perfect fluid is inequivalent to the coupling to another massless 
scalar field % in a quartic way g'^cjy^x^- If we study such a system: 



inG 



8ttG 



Px\ 



(12) 



By considering = /2 + g'^^P'x^ the structure of the corresponding interacting term is different from the 
one in Eq. 

As a final comment of this section, it is useful to compare this setting with the T = case, which is very similar to 
the single inflaton case since pp oc a"^^^^^''^ is rapidly washed out for wp > 0. The dissipative term F slows down 
the inflaton during the slow-roll regime and damps its amplitude during the oscillatory stage. As an example, we 
show the left panel of Fig. [l]the evolution of a massive inflaton - ¥{4)) = m^0^/2 - with F = 0.5m. See the caption 
for the evolution for other relevant background quantities. 



III. SCALAR PERTURBATIONS IN THE LONGITUDINAL GAUGE 

The scalar perturbations around a flat Robertson- Walker metric are: 

ds^ = -(1 + 2a)df ~ af3^^dtdx' + c? [5^^ (1 - 20) -f 27,^] dx'^dx^ , (13) 
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where the symbol ,i denotes the derivative with respect to the spatial coordinates. 
If we work in the longitudinal gauge, the perturbed metric is: 



(14) 



where the two longitudinal perturbations are the same since neither the scalar field nor the perfect fluid have linear 
anisotropic terms in the pressure. The energy-momentum tensor for the fluid is: 



Tf_,„ ^ {p + p)ui_,Uy +pg^^ , u^M^ = -l 

and its perturbations are given by: 

= SpF , STpi = ipF ,1 STpj = SpF Sj = wpSpF 5] , 

where we have implicitly introduced the potential ipF for the spatial velocity of the fluid Ui 
The perturbed total energy and pressure are, respectively: 

Sptot = (/> (50 - + V^S(I) + 6pF , 



(15) 
(16) 

(17) 



Sptot = (t)S4> ~ - V^64> + lofSpf ■ 



The Einstein constraints are 



fc2 

3i7$ + 3ij2$ + _$ = -AnGSpt. 



(18) 



(19) 



i7$ + $ = AttG (t)5<t) - i/^F 



and the conservation of the energy momentum tensor for the fluid leads to: 

6pF + SH{1+u;f)5pf 5-^/'i? - 3(1 + w_F)/9i^$ = Xf'^ + SXf 



= 2r 



IpF + SHipF + $(1 + ijJF)pF + UJFSpf 



with 



SXf = 2T^5(t) - 2r(/)2$ + 02 = -SX^ 
The equation of motion for field fluctuations is: 



{3H + r) , 



4 0$-2K6$-r04'-0(5r 



The curvature perturbation in the longitudinal gauge C is defined as: 

H\ J 3H{l + cutot) 

Following J22'] we define the entropy perturbation between the scalar field and the fluid: 



S. 



4>F 



_ 3H^Fl,pPF f Sp^ Sp 



1F14>PF 



Spf 



7totPtot \ P4, Pf ) 7totPtot VlFpF - r^4,P4, 70^0(1 + 7") 
and the intrinsic entropy perturbation of the scalar field: 



^Hj4,cl / Sp^ 



l-ci 



clSp4> 



p4l-4)(l + r) 



(20) 



(21) 
(22) 

(23) 
(24) 

(25) 

(26) 

(27) 
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We conclude this section with the equation of evolution for C: 



c 



3^ — $ + 47rG(5pnad 



(28) 



where (5pnad = ^Ptot — c^Jptot is the non-adiabatic pressure which depends on the intrinsic and relative entropy 
perturbations, as we can see making use of Eq. ([5]): 

^'''^^^ =(l-4)5, + (c.p-4)(l-^r)5^,. (29) 



P0(l + r) ^ 0/ ■ V - ^^p^ 

It is simple to see from Egs. ipS)) and that in the case of pure radiation (0 = = = and LOp = = 1/3) 
curvature perturbations are conserved on large scales. 

Large Scale Curvature Solution during Slow-Roll 

It is useful to give the analytic approximation of the adiabatic solution for large-scale cosmological perturbation 
during slow-roll. By using Eqs. (|9|10llip . the approximate solutions at leading order in the slow- roll parameters are: 

$ ~ —— 0(1-1- r),5(/) 
tl 

ipF — —r(l)S4> , 



oc 



|, (30) 



with j3 ~ f/(3iJ^). These solutions have been obtained neglecting the highest time derivatives and the terms oc in 
Eqs. (|20l21l22l24p and are valid if F = r(0). It is interesting to use the above relations and write the leading term 
(in slow- roll parameter) of curvature perturbations: 

2 $ $ 

C-:tt- --^^50, 31 

3 1-1- cjtot e 

which is formally just the single field expression; T is however present in the dynamics for the inflaton, which differs 
from single field case with the same V{(j)) and F = 0. The expression differ from the multi-field inflationary case in 
which all the scalar fields in slow-roll contribute (with the same formal weight) to curvature perturbations. We end 
on noting that C, is constant in time for the last expression of Eq. ((30)) at leading order in slow-roll parameters (i.e. 

C ~ ffCO(e))- 



IV. NUMERICAL EVOLUTION OF COSMOLOGICAL PERTURBATIONS 

We will consider the numerical evolution of linear cosmological perturbations from inflation through reheating, 
following Ref. [5|] which evolved the Mukhanov variables related to two scalar fields in a fully constrained evolution. 

In this case we need to obtain the evolution of the Mukhanov variable associated to the scalar field coupled with 
Lukash variable associated to the gradient velocity of the perfect fluid. Given a fully general metric perturbation as 
in Eq. p^ . gauge invariant variables for the scalar field and the fluid are defined as: 

= <50 + 4't ^^F-{p+ P)i^n ■ (32) 

In the UCG, where the spatial metric is unperturbed, 



ds^ = -(1 + 2a)dt^ - ajS^idtdx^ + a'^Sijdx'^dx^ , 



(33) 
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the gauge invariant variables coincide with the scalar field and fluid velocity potential fluctations. It is therefore easier 
to obtain the evolution for these gauge invariant variables in such a gauge. From here on we restrict ourselves to this 
gauge and we shall write fluctuations of the scalar field and of the fiuid, although different, in the same way as the 
previous section. 

The equations corresponding to (|21I22I24[) are: 



J^2 ^2 
5pF + '^H{1 + LOf) SpF jV'F ^ (1 + ^f)/3_F7P/9 — XFCt + SXF 



2a 



= 2r (P S4> - Tcp^ a + ST < 



(34) 



ipF + 3HipF + ct{l + ujf)pf + ujfSpf = ~TipS4>. 



(35) 



{3H + r) ( 



fc2 



Scf) = (f)[a + —(3 ) - 2V^a - (j)ST - aVcj) . 



(36) 



If S(f> is used to quantize the linear scalar field fiuctuation in presence of gravity [20] , the variable to quantize a fluid 
[2l| . as also explained in [2^, is: 



Qi 



ijjf 



^/pf{1 + wf) 



^ij Qj — 



where G has the following components 

G^^ — 

G^F = 

Gfa — 



47rG(w| - 1)pf0 
y/pf(1 + wf)wfH 
[-47rG'(w| - l)pF + {2wF + l)HT] 



y/pF{l+WF)H 



Gf f — r — 

PF 



(37) 



On restricting to the case in which F is constant in time, the coupled second order differential equation for {Q^ , Qf) 
{54> , Qf) are: 



(38) 

(39) 
(40) 

(41) 

(42) 
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and where f2 is given by 



16.^G^i^('p.(^^^il±ll! + 20^) (43) 



3 (wi. - 1) r i \ 



-327r2G2^ + 27rG^^^i ^77— - 47rGr^ I (44) 

H'^ Wf ti pf ti I 

10 , ^i , (3i(;Fi? - r) r,/. 



f^F0 = Vpf(1 + wf) "4(1 + WF)vrG-^ - 127rG(l + wf 



H {1 + wf)pf (l + Wi^)Pi 



+327r2G2(l + u;^)^ + 16tt'G'{1 + wf)^ + ^^G , ^ | (45) 



J^FF = WF^ + 7g'(l-^^F) + 6^Gpf(l + u;F)(l + 3wf)- '^^ {I + wfY 

—1 — 1- 67rG0 (w_F — 1) — h - — (—21 + 6wFti) 

Pf H 2pF 

^^PF0 (l + ii'F) (46) 

Several cross-checks on the novel set of equations can be made: 

1) For Pf = F = the equation for agrees with the Mukhanov equation for single field inflation [20l |. 

2) When the two components are decoupled (F = 0) and the fluid is stiff matter (wf = 1) the above equations 
agree with those for two scalar fields Q when the second one x is decoupled and massless (i.e. equivalent to stiff 
matter) reminding that = —Qf- 

We now evolve the system in Eq. ([38|) starting from wavelength well inside the Hubble radius fc >> aH with initial 
conditions: 

Qf - -777^-TT7I ■ (47) 



~ a(2fc)i/2 ' a(4fc2u;f.)i/4 
The gauge-invariant curvature perturbation C, in UCG gauge is written as: 



C = i^^%^^^^^^^^C. + CF. (48) 
4>^ + Pf[1 + Wf) 

Evolution of curvature perturbations are displayed in Fig. [2l Perfect fluid fluctuations are amplified similarly to 
inflaton fluctuations when F 0, leading to a mixture of curvature and isocurvature fluctuations during slow-roll. 
Isocurvature perturbations can then change the amplitude curvature perturbations during and after the inflaton decay. 

The uncoupled case 

As is clear from the numerics, the perfect fluid does not contribute to curvature perturbations in the uncoupled case 
F = 0. We have already observed that for wf > 0, the homogeneous perfect fluid energy density is rapidly washed 
out for F = 0. If we neglect the terms which decay quasi-exponentially on time, the equation for fluid perturbations 
becomes: 



Qf + 3HQf 



fc2 g . 
WF^ + -H^{1 - wl) -f QTlG(j)'^(wF - 1) 

4 



. (49) 



For radiation, at leading order in a nearly de Sitter exponent, we have an effective mass ~ 2H^ which makes the fluid 
fluctuations similar to a massless conformally coupled fluid. The initial vacuum fluid fluctuations in Eq. (|47p are 
almost left unchanged during the streching to large scale, leading to a spectrum ^^(fc) = k'^\Q f]"^ I {2ii^) ^ fc^, which 
is too blue to affect the nearly scale-invariant spectrum of inflaton fluctuations in the observable range. 
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FIG. 2: Typical evolution of curvature perturbations (total as solid line, short-dashed and long-dashed for the field ("^ and fluid 
contribution, respectively) from inflation to radiation as function of an adimensional cosmic time mt. On the left panel, 
the evolution for T/m = 0.5 for which the curvature perturbations are given by the inflaton (solid and short-dashed lines are 
super-imposed) and there is negligible post-inflationary change. On the middle and right panel, the evolution for T/m — 6 and 
r/m = 8. 

V. DISCUSSIONS AND CONCLUSIONS 

We have studied the coupling of the inflaton to a perfect fluid through a friction term F during and after inflation. 
Although this type of coupling has been used in the early days to describe the decay of the inflaton after the accelerated 
era, it has also been used in the regime T ^ H for warm inflation. 

By considering radiation as the perfect fluid and F constant in time for simplicity, we have shown how this coupling 
freezes the perfect fluid during the slow-roll evolution and large-scale fluid fluctuations are amplifled together with 
the scalar field ones, leading to a mixture of curvature and isocurvature fluctuations during the slow-roll regime. In 
the stage of the inflaton decay, large scale curvature and isocurvature perturbations are not weakly coupled anymore 
and there is a transfer of the latter leading to an amplitude for the former which is different from the one computed 
during slow-roll. After the decay is completed, a purely adiabatic curvature perturbation is left. We believe that the 
qualitative behaviour found here is generic and does not occur only for the simplest parameters studied in this paper 
- radiation as the perfect fluid and F constant in time. Note that a variation of large scale curvature perturbations 
was also found in the late decay of a massive curvaton [23 |. 

Our results for a coupling to a perfect fluid add to the changes of curvature perturbations during preheating found 
for other couplings to scalar fields, as studied in As already said, this coupling of the inflaton easily generates 
non-negligible isocurvature fluctuations during inflation, which is one of the key requirements in order to have a 
change of curvature perturbations in the post-inflationary stage 0. Since we see that the amount of isocurvature 
perturbations increases with F in the cases studied here, inflationary models with sizeable r have a larger variation of 
curvature perturbations in the post-inflationary era. As shown explicitly here, a successfuU conversion of isocurvature 
into curvature perturbations at reheating may not need a varying decay rate. 
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